High-precision measurement of the proton elastic form 
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Abstract 

We report a new, high-precision measurement of the proton elastic form 
factor ratio UpGe/Gm for the four-momentum transfer squared = 0.3- 
0.7 (GeV/c)^. The measurement was performed at Jefferson Lab (JLab) in 
Hah A using recoil polarimetry. With a total uncertainty of approximately 1%, 
the new data clearly show that the deviation of the ratio ^pGe/Gm from unity 
observed in previous polarization measurements at high continues down to 
the lowest value of this measurement. The updated global fit that includes 
the new results yields an electric (magnetic) form factor roughly 2% smaller (1% 
larger) than the previous global fit in this range. We obtain new extractions 
of the proton electric and magnetic radii, which are {r'^y^'^ = 0.875 ± 0.010 fm 
and {r\,.j)^^^ = 0.867 ± 0.020 fm. The charge radius is consistent with other re- 
cent extractions based on the electron-proton interaction, including the atomic 
hydrogen Lamb shift measurements, which suggests a missing correction in the 
comparison of measurements of the proton charge radius using electron probes 
and the recent extraction from the muonic hydrogen Lamb shift. 



The nucleon electromagnetic form factors are fundamental quantities which 
relate to the charge and magnetization distributions within the nucleon. For 
over 40 years, the form factors have been studied extensively by Rosenbluth 
separations of the unpolarized electron-proton scattering cross section. Recent 
advances in the technology of intense polarized beams, polarized targets and 
polarimetry ushered in a new generation of experiments that measure double 
polarization observables [H, [3j S] ■ Although the proton electric to magnetic form 
factor ratio R = UpGe/Gm determined by unpolarized measurements showed 
minimal dependence up to « 6 (GeV/c)^, experiments at JLab with high- 
quality polarized electron beams measuring recoil polarization 0, 0, Q revealed 
that the ratio UpGe/Gm drops almost linearly with Q'^ for Q^>1 (GeV/c)^. 
These findings led to an explosion of experimental and theoretical efforts to 
understand the proton form factors The difference between the po- 

larization and cross section measurements is now believed to be the result of 
two-photon exchange (TPE) contributions 0, [13, [Hi , which have little im- 
pact on the polarization measurements but significantly affect the Rosenbluth 
extractions of at large Q^. 

While measurements at large momentum transfer have provided information 
on the fine details of the proton structure and relate to the quark orbital angu- 
lar momentum, the low form factor behavior is sensitive to the long-range 
structure which is believed to be dominated by the "pion cloud" . High-precision 
measurements at low were motivated by a recent semi-phenomenological 
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Table 1: Kinematic settings for the experiment. 9e (8p) is the scattered electron (proton) angle 
in the lab frame, Pp is the proton central momentum, e is the virtual photon polarization. 



Q2 




Op 


Pp 


e 


Analyzer 


}eV/c)2 


(deg) 


(deg) 


(GeV/c) 




thickness 


0.298 


28.5 


60.0 


0.565 


0.88 


2.25" 


0.346 


31.3 


57.5 


0.616 


0.85 


2.25" 


0.402 


34.2 


55.0 


0.668 


0.82 


3.75" 


0.449 


36.7 


53.0 


0.710 


0.80 


3.75" 


0.494 


39.2 


51.0 


0.752 


0.78 


3.75" 


0.547 


41.9 


49.0 


0.794 


0.75 


3.75" 


0.599 


44.6 


47.0 


0.836 


0.72 


3.75" 


0.695 


49.8 


43.5 


0.913 


0.66 


3.75" 



fit [12| . which suggested that structure might be present in all four nuclcon 
form factors at ~ 0.3 (GeV/c)^. Later polarization measurements from 
MIT-Bates [Hi and JLab [3 probed this region with a precision of ^2%, but 
found no indication of such structure in the ratio ^pGe/Gm- 

This Letter reports on a new, high-precision polarization transfer measure- 
ment (JLab E08-007) of the proton form factor ratio UpGe/Gm at values 
between 0.3 and 0.7 (GeV/c)^ . In the one-photon exchange (Born) formalism, 
the ratio of the transferred transverse to longitudinal polarizations is related to 
the proton form factors 0: 

Ge Ee + E' 9e.Pt ,1 X 

where Mp and fip are the proton mass and magnetic moment, Ee [E'^) is the 
incident (scattered) electron energy, 0^ is the electron scattering angle and Pj 
{Pi) is the transverse (longitudinal) component of the polarization transfer. 



The experiment was performed at Jefferson Lab in Hall A [15[. A 1.2 GeV 
polarized electron beam of between 4 and 15 /iA was incident on a 6 cm thick 
liquid hydrogen target. The beam helicity was changed at 30 Hz, with a quartet 

structure selected pseudorandomly between (H h) and ( — h +— ) for each 

set of four helicity states. The beam polarization was near 83%, as measured 
by the M0ller polarimeter in the Hall [l^. The recoil proton was detected by 
the left High Resolution Spectrometer (LHRS) in coincidence with the elasti- 
cally scattered electron detected in a large acceptance spectrometer ( "BigBite" ) . 
The transferred proton polarization was measured by a focal plane polarimeter 
(FPP) in the LHRS [l5|. The trigger was a coincidence between detection of 
a charged particle in the HRS and a signal in the BigBite calorimeter used to 
select energetic electrons. Since elastic events can be well identified using the 
proton kinematics, only information from the BigBite calorimeter was used in 
the analysis; the tracking detectors were not turned on for the experiment. The 
kinematic settings are given in Table [1] 

The FPP measured the azimuthal asymmetry due to the spin-orbit coupling 
in proton scattering from carbon nuclei. The proton's transferred polarization 
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was extracted by a weighted-sum method [16| based on a COSY [17| model 
of the spin precession in the spectrometer. Due to the fast reversal of beam 
helicity, the helicity-independent corrections (acceptance, detector efficiency, 
target density, etc.) cancel. At each value, measurements were taken at the 
nominal elastic kinematics and with the HRS spectrometer momentum shifted 
by +2% and -2% to allow for additional systematic studies of the spin transport 
model. The kinematic factors in Eq. [1] were determined from the beam energy 
Ef, and the proton scattering angle Qp. 

A series of cuts were applied to cleanly select the elastic events and minimize 
the systematic uncertainties. In addition to the standard HRS cuts (tracking, 
acceptance, timing, etc.), cuts on the FPP variables were applied to select the 
events with reliable second-scattering reconstruction [l^l- For elastic event se- 
lection, we use the correlation between the proton angle and momentum. Af- 
ter applying a tight cut (l.Ta) on the elastic peak and a target vertex cut, 
the contamination from the aluminum endcaps of the target and neutral pion 
photoproduction is less than 0.1% in the extracted ratio /ipG^/Gi/, which is 
negligible compared to other systematic uncertainties. 

The primary systematic uncertainty in the measurement comes from un- 
certainty in the spin precession of the protons in the spectrometer magnetic 
fields. The uncertainty arising from imperfect knowledge of the spectrometer 
fields was estimated by examining the change in the results when key param- 
eters of the spectrometer optics model used for spin precession were altered. 
The most sensitive parameter is the dipole central bending angle 6o, which is 
nearly 45°. A conservative 5.5 mrad uncertainty of this parameter is quoted 
from previous analysis yielding a 0.1-0.6% change in the ratio for different 
kinematics. In addition, any error in the particle trajectory reconstruction will 
lead to incorrect proton kinematics at the target which are used as input to the 
spin precession calculation. The largest effect comes from offsets in the proton 
angle in the scattering plane, estimated to be ~1 mrad based on the consis- 
tency of the particle kinematics with the elastic constraint when accounting for 
the uncertainties of other related parameters (beam energy, proton momentum, 
spectrometer angle). The associated uncertainty on the extracted form factor 
ratio is determined to be 0.6-0.9% by manually shifting the proton trajectories 
in the spin transport [l^. As a final check, the extracted value for [i^GeIGm 
was examined as a function of each of the reconstructed proton target variables 
and for each of the three momentum settings taken for each setting. All of 
these results were consistent within the statistics of the measurements. 

The experimental results are summarized in Table [5] The average FPP 
analyzing power (Ac) and efficiency e arc in good agreement with previous mea- 
surements [l^. When a transversely polarized proton interacts with the FPP, 
the azimuthal asymmetry induced by the secondary scattering is proportional to 
the proton polarization and the analyzing power, which depends on the proton 
momentum and the angle of the secondary scattering. The average analyzing 
power extracted in the analysis is the event-weighted average over the 5-25° 
range of angles used in the analysis of the experiment. The FPP figure of merit 
(FOM) is an integral of over the selected Q range. No correction was ap- 
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Table 2: Experimental ratio fipG e /G m it stat. it syst. along with the average FPP analyzing 
power {Ac), efficiency e and the figure of merit {eA^) for secondary scattering angles between 



5° and 25° 



Q2 


(Ac) 




FOM 




(GeV/c)2 




[%] 


[%] 




0.298 


0.219 


5.30 


0.25 


0.9272±0.0114±0.0071 


0.346 


0.394 


3.67 


0.57 


0.9433±0.0088±0.0093 


0.402 


0.466 


4.36 


0.95 


0.9318±0.0066±0.0076 


0.449 


0.488 


4.09 


0.97 


0.9314±0.0060±0.0073 


0.494 


0.466 


3.81 


0.83 


0.9286±0.0054±0.0076 


0.547 


0.430 


4.34 


0.81 


0.9274±0.0055±0.0071 


0.599 


0.392 


4.41 


0.68 


0.9084±0.0053±0.0104 


0.695 


0.334 


4.74 


0.53 


0.9122±0.0045±0.0107 



plied for two-photon exchange, which is less than a 0.1% correction on Pt/Pi 
(or Ge/Gm) for these kinematics [§]. 

With the unprecedented statistical precision, we were able observe a small 
cut dependence in the analysis when loose cuts (larger than 2cr) were placed 
on the elastic peak, and so used tighter cuts (1.7(t) than previous experiments 
where were statistics limited. This also led to a detailed reanalysis of a previous 
experiment IJ], which found that their sensitivity to these effects were small 
but yielded an improved isolation and correction for events involving scattering 
from the cryotarget endcaps that led to a slight decrease in the extracted form 
factor ratio [201 . 




Figure 1: Low polarization measurements of ^pGe/Gm- The solid red circles are results 
from E08-007 and the hollow symbols are from previous measurements \T3l. |2C|. |21| with 
precision of better than 3%. The grey crosses show various lower precision data J22,[23,|2J, 2^. 
The solid curves are selected theoretical models of the proton form factors I26ll27l . |2S|, |29|, 30| . 
The dashed curves are the Friedrich and Walcher fit 11211 . and the updated global fit (see text). 
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Figure [T] presents the new results together with previous polarization mea- 
surements, including the updated results from Ron, et al [l^. The inner error 
bars are the statistical uncertainty and the outer error bars are the statistical 
and systematic uncertainties added in quadrature. The new results achieved a 
total uncertainty of ^--^1%, usually dominated by the systematic uncertainties. 
The data fall slowly across the covered range and strongly deviate from unity. 
The high-precision point at = 0.8 (GcV/c)^ [2l[ is also consistent with the 
trend of our new results. Compared to the Bates measurement these new 
ratio results are lower by 2-3cr, indicating some systematic discrepancy. 

Two implications of the new data can be directly seen. Some previous mea- 
surements suggested that the ratio was flat at fipGE/GM ~ 1 until « 
0.2 (GeV/c)^ and then began to fall, which implied identical low-Q^ behav- 
ior of the form factors and thus equal charge and magnetization radii for the 
proton. We observe that ^jLpGe/Gm is significantly below one over the entire 

range of our data, with no indication that the ratio is approaching unity no- 
ticeably above = 0. The global fit (described later) suggests that the ratio 
begins to fall starting at or very near = 0, implying significant differences in 
the large scale spatial distribution of charge and magnetization in the proton. 
Second, with this high precision data, we see no indication of any structure in 
the ratio ^pGe/Gm over the covered region. 

The curves in Fig. [T] show a selection of modern fits and models in the 
literature: a light-front cloudy-bag model (Miller [1^), a point-form chiral con- 



stituent quark model (Boffi et al. 27[), a Lorentz covariant constituent quark 



model (Faessler et al. 12811 ) . and two representative vector- meson dominance 
models (Belushkin et al. |29l| and de Melo et al. js^l ) • None of the existing mod- 
els precisely match the new results, but we note that despite utilizing different 
theoretical frameworks, both the calculation of Miller and of de Melo et al. are 
in qualitatively good agreement with our results. As both of these consider the 
role of non-valence contributions, e.g. nucleon plus pion contributions in addi- 
tion to the lowest 3-quark Fock state, this provides additional support for the 
important role of the pion cloud at low Q^, even though no structure is seen in 
this ratio (or in the global fit for Ge and Gm)- 

We have performed an updated global fit of the form factors, following the 



procedure of Ref. [3]|. In addition to the data sets included in that anal- 



ysis, we add our new results along with other recent polarization measure- 



ments [y, l20|, l2l| , but have not included the more recent cross section measure- 



ments from Mainz [32| , as these data do not include sufficient information on the 
uncertainties to be included in a consistent fashion. Figure [2] presents Ge and 
Gm normalized by the standard dipole form, Gd = + Q^/0.71(GcV/c)^]^, 
for the previous global fit and our updated fit, along with the Rosenbluth sep- 
aration value for the cross section analysis of Ref. [3l| without the inclusion 
of TPE corrections. The difference between the cross section analysis (hollow 
circles) and the global fit (dotted curve), both from Ref. [3l|, is due to a com- 
bination of the TPE corrections and the inclusion of the limited polarization 
data available at the time. The difference between the dotted curve and our 
updated fit (solid red curve) comes from the inclusion of the new polarization 
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data (Refs. [H [2^ and this work). For 0.2<Q2<i.o (GeV/c)^ the updated 
global fit shows a 2% decrease in Gb and a 1% increase in Gm compared to the 
earlier fit (sj . The inclusion of these precise polarization data also improves the 
determination of the relative normalization of the different cross section data 
sets included in the fit, which has a significant impact on the uncertainty in the 
extraction of the proton charge radius . 




[GeV/cf 

Figure 2: Global fits to Gb, Gm and ratio (IpGe/Gm from ref. |31|| (dotted line) and this 
work (solid red line). The blue lines are the best fit and uncertainty band from the Mainz 
cross section data [32 |. The solid circles are our new results for ^IpGe/Gm and the hollow 
circles show the global cross section analysis of Ref. [31II . which do not include TPE TPE 
corrections. The green long-dashed line shows the first term of our low-Q^ fit used to extract 
the radius. 



We also show the fit to recent Mainz cross section measurements [34I in 
Fig. [5] The Mainz fit is in reasonable agreement with our direct extraction of 
Mp^b/Gmj but their results for Ge and especially Gm are somewhat incon- 
sistent with our global fit (and systematically above all previous extractions of 
Gm)- However, it is difficult to make a detailed comparison of the results due 
to differences in TPE corrections. The hollow points do not include any TPE 
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correction, the global fits include the full TPE corrections of Ref. Q to the cross 
sections, while the Mainz extraction [slj applies only the = limit of the 



Coulomb distortion correction [3J] . This yields an overestimate of the correction 



at all finite values and neglects the dependence of the correction which 



will change the extracted magnetic form factor and magnetic radius [35l l36l . |2(J | , 
but has a much smaller impact on the very low behavior of G^; and the ex- 
tracted charge radius. A full examination of the consistency between the Mainz 
result and other cross section measurements would require careful examination 
of the TPE corrections applied to the different data sets and is beyond the scope 
of the present work. 

The proton electric and magnetic RMS radii can be determined from the 
form factor slope at = 0: 



We do not use the previously described global fit to extract the radius because 
such fits are dominated by high-Q^ data which have essentially no sensitivity to 
the radius. To extract the radius, we repeat the global fit described above using 
a continued fraction [s^ parameterization with fewer parameters, including only 
cross section and polarization data below (GeV/c)^. The radius values 

we obtain are: 

(r|)i/2 ^ 0.875 ± O.OOScxp ± 0.006fit fm (3) 
{rliV^^ = 0.867 ± 0.009cxp ± O.OlSfit fm, (4) 

where the quoted values are the averages from fits with 3, 4, and 5 parameters, 
the first uncertainty is the result of the experimental uncertainties (statistical, 
systematic and normalization uncertainties added in quadrature) and the second 
is the model dependence estimated by taking the scatter of fits when we vary 
the functional form, number of fitting parameters, and cutoff. The slope 
corresponding to first term of the low-Q^ analysis is shown in Fig. [2j For Ge, 
the global fit and low-Q^ fit arc in good agreement as 0, although by 

= 0.1, the higher-order terms are beginning to yield a deviation from the 
linear term. For Gm, even the slope at = is slightly different. Because the 
cross section becomes more sensitive to Gm as increases, fits that include 
high data will tend to overemphasize agreement at high values at the cost 
of the low behavior. This is why the global fits are generally not well suited 
to extraction of the radii, especially when fitting the magnetic form factor. 

While the charge and magnetic radii yield consistent values, as one would 
expect in the non-relativistic limit where the quarks carry both the charge and 
magnetization, the fit yields a different detailed behavior of the form factors 
even at low . This implies a different shape for the charge and magnetization 
distributions, even though the radii are consistent. If a fixed functional form is 
chosen for both the charge and magnetic form factors at low Q^, the fits yield 
different radii. 



8 



Our magnetic radius is significantly larger than the Mainz value [s^I : = 
0.777 ± 0.017 fm. As noted previously, the Coulomb correction applied in that 
analysis does not include any dependence, which will likely lead to an er- 
ror in the extracted magnetic radius. An analysis of their data with a modern 
correction for Coulomb distortion (or full two-photon exchange) is required to 
determine if this represents a true discrepancy between the Mainz result and 
all other extractions. An initial evaluation [ssl, [s^, using a calculation for the 
TPE corrections valid up to « 1.0 GeV^, yields an increase in the magnetic 
radius of 0.026 fm, accounting for roughly one-third of the discrepancy. 



Sick 






Bernauer et al. ; 
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— •- 




CODATA 


■— 








Poll! et al. * 

' • 







0.82 0.84 0.86 0.88 0.90 
Proton charge radius [fm] 



Figure 3: The proton RMS charge radius from a previous ep scattering analysis (Sick l33t ). the 
Mainz low-Q^ measurement (Bernauer et al. [s^ l ) a nd this work compared to the CODATA [37j 
and muonie hydrogen spectroscopy (Pohl et al. Km ). The red dashed hncs show the combined 
results from CODATA, Bernauer et al. and this work while the black dotted lines show the 
Pohl et al. uncertainty. 



Several extractions of the proton charge radius are summarized in Fig. [3l 
Based on totally different data sets, the charge radius extracted here is in 
excellent agreement with the Mainz extraction 0.879 ± 0.008 fm [s^l and the 



CODATA value 0.8768 ± 0.0069 fm [37[, which is mainly inferred from atomic 



hydrogen spectroscopy. The combined result of the three measurements based 
on electron-proton interactions gives: 



{rl)\/p^ = 0.8772 ± 0.0046 fm. (5) 
The Sick result is shown for comparison but not included in the average because 



our updated extraction includes all of the data included in Ref. [33[. A recent 
measurement of the Lamb shift in muonie hydrogen (ssj yields a proton charge 
radius of 0.8418 ± 0.0007 fm. which is 7.7(t from the combined result of the 
ep measurements. Barring an error in the extraction of the muonie hydrogen 
result or a common offset affecting both electron scattering and atomic hydro- 
gen measurements, there appears to be a true difference between electron-based 
and muon-based measurements. This would seem to disfavor explanations that 
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would lower the CODATA result, e.g. based on a modified value of the Ryd- 
berg constant, and suggests a missing correction that differs for electron and 
muon probes. One mechanism that could cause such a difference was recently 
proposed [s^. However, there is still significant activity aimed at determining 
whether the discrepancy is the result of a missing element in the QED correc- 
tions in the bound fip system (40| or if a more exotic explanation is required. 

In conclusion, we present a new set of high precision form factor ratios 
extracted using the recoil polarization technique. These data provide improved 
sensitivity to Gm at low while minimizing two-photon exchange corrections 
which are of particular importance in separating Ge and Gm- Our results 
have been combined with previous cross section and polarization measurements 
to provide an improved global analysis of the form factors. The new data 
shift the extracted values of Ge and Gm at low Q^, which will impact the 
proton structure corrections to atomic energy levels in hydrogen and muonic 
hydrogen as well as the extraction of the strange quar k contribution to the form 
factors probed in parity- violation ep scattering |4ll . |4^ . While the updated 
form factors yield only a small modification to the expected parity-violating 
asymmetry, typically less than half of the assumed uncertainty, it will have a 
correlated effect on the analysis of all such measurements. Finally, we present 
an improved extraction of the proton charge and magnetization radii, and find 
consistency between electron-based probes of the proton radius, which disagree 
with the muonic hydrogen result. 
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